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Abstract 

Diagnosing core of supernova requires favor-dependent reconstruction of three species of neutrino 
spectra, i> e , v e and v x (a collective notation for Vn, i'u, v T , and P T ). We point out that, assuming 
the information available, CPT symmetry can be tested with supernova neutrinos. We classify all 
possible level crossing patterns of neutrinos and antineutrinos into six cases and show that half 
of them contains only the CPT violating mass and mixing patterns. We discuss how additional 
informations from terrestrial experiments help identifying CPT violation by narrowing down the 
possible flux patterns. Although the method may not be good at precision test, it is particularly 
suited to uncover gross violation of CPT such as different mass patterns of neutrinos and antineu- 
trinos. The power of the method is due to the nature of level crossing in supernova which results in 
the sensitivity to neutrino mass hierarchy and to the unique characteristics of in situ preparation 
of both v and v beams. Implications of our discussion to the conventional analyses with CPT 
invariance are also briefly mentioned. 

PACS numbers: 11.30.Er,97.60.Bw,14.60.Pq 



*Electronic address: E-mail: minakata@phys.metro-u.ac.jp 
^Electronic address: E-mail: uchinami@phys.metro-u.ac.jp 

Typeset by REVTeX 1 



I. INTRODUCTION 



CPT is one of the most fundamental symmetries in relativistic quantum field theory 
by which masses and flavor mixing angles are constrained to be identical for particles and 
their antiparticles. Because of the fundamental nature of the symmetry, it is important to 
test the CPT invariance, and there has been continuing efforts mainly in kaon physics 0. 
Needless to say, the effort of testing CPT symmetry should be extended to the lepton sector. 

Recently, there arose some interests in possible violation of CPT symmetry in the lepton 
sector, partly motivated by a possible interpretation of the LSND result .3]. A hypothesis of 
different mass and mixing patterns of neutrino and antineutrino sectors, as first suggested 
by Murayama and Yanagida [4] , could be flexible enough to accommodate the LSND data in 
the three-neutrino framework, while not sacrificing the success of describing the atmospheric, 
the solar, the reactor, and the accelerator data [a. El 0, HI • While the proposal was followed 
by a series of papers 0,0, 11 , it was shown by Gonzalez-Garcia, Maltoni and Schwetz 12l |. 



in an extensive statistical analysis of all the data including KamLAND [j], that the CPT 
violating hypothesis is not in good shape. Interestingly, the best fit point of all data except 
for LSND is CPT symmetric, and the mixing parameter region favored by LSND is more 
than 3<T away from the region favored by all but LSND data. 

In this paper, we explore a possibility of testing CPT symmetry with supernova neutrinos. 
Independent of the success or the failure of the CPT violating scenario for the LSND data, 
it is important to test CPT as a fundamental symmetry in nature. Given the fact that 
neutrinos has brought us several surprises, there exists an even more intriguing (albeit not 
likely) possibility to discover CPT violation by future neutrino experiments. Supernova 
neutrinos are advantageous to examine neutrino and antineutrino properties simultaneously 
and consistently because the beam is composed not only of z/ e , z/ M , and v T but also their 
antiparticles. We examine the possibility of using neutrinos from supernova to identify 
CPT violation assuming the resolving power of flavor-dependent neutrino fluxes in future 
observation of galactic supernovae. 

So far, constraints on CPT violation of mixing parameters in the lepton sector have 
been derived by Super-Kamiokande (SK) group |l3| and in jl2|. Possible ways of testing 
CPT symmetry has been discussed by using solar and reactor neutrinos 3 Hl| > and neutrino 
factory [16]. We restrict ourselves, as these preceding works do, to the framework of possible 
CPT violation in masses and flavor mixing of neutrinos, assuming that neutrino interactions 
conserve CPT. Then, the natural question is how supernova method can be competitive 
to these more "traditional" methods for testing CPT. It is the right question because it is 
unlikely that such rare event as supernova can be used for a precision test of CPT symmetry. 
Despite the reasonable skepticism, we will show in this paper that supernova neutrino can 
be a powerful tool for uncovering gross violation of CPT symmetry. (See Sec. [H]for further 
comments.) 

We rely on the reference for the formulas of neutrino flavor conversion in supernova in 
the three-flavor framework. We should note that our general CPT non- invariant treatment, 
of course, includes the case of CPT invariance. Therefore, the reader can use part of the 
formulas given in this paper as a compact recollection of those in but by now without 
ambiguities due to the solar neutrino solutions. 

In Sec. HH we discuss the question of why and how supernova neutrinos are useful to test 
CPT. In Sec. IIII[ we review the basic properties of supernova neutrinos and the approxima- 
tions involved in our treatment. In Sec. II VI after recollecting compact formulas for neutrino 
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flavor conversion in supernova, we present a complete classification of spectral patterns of su- 
pernova neutrinos that are possible in a general CPT violating ansatz. In Sec. we discuss 
how the allowed flux patterns of supernova neutrinos reduce as additional input of 6*13 and 
neutrino mass hierarchy are added. In Sec. IV11 we give a comparative study of characteristic 
features of spectra of three effective neutrino species predicted in each classified patterns 
of neutrino flavor transformation in supernova. In Sec. IV11[ we discuss at a qualitative (or 
semi-quantitative) level to what extent the possible different neutrino flux patterns can be 
discriminated observationally. In Sec. IVIIIl we give the concluding remarks. 

II. WHY AND HOW ARE SUPERNOVA NEUTRINOS USEFUL TO TEST CPT? 

To answer the question of why and how supernova neutrinos are useful to test CPT we 
need to specify the question; Namely, which aspects of CPT symmetry do we want to test, 
or which features of CPT violation do we try to uncover? 

First of all, lacking well defined models of CPT violation, we cannot test it in neutrino 
interactions by using supernova neutrinos. If the interactions are different from what we 
know the neutrino properties inside the core must be recomputed with new interactions to 
define CPT non-invariant features of supernova neutrinos. We do not have the recipe to 
carry it out. Therefore, we restrict the type of CPT violation to test to the ones signaled by 
difference between masses and mixing parameters of neutrinos and antineutrinos, assuming 
that their interactions are described by the standard model. 

How can CPT violation be actually signaled by supernova neutrinos? As we will show in 
the subsequent sections, possible difference in mass patterns and mixing angles of neutrinos 
and antineutrinos results in several different spectral patterns of three species of neutrinos, 
v ei v ei and u x , where the last is a collective notation for u^, v T , and v T . (See Sec. IIIII ) 
In this paper, we rely on the assumption that flavor-dependent reconstruction of supernova 
neutrino fluxes will be done at the time of next supernova, so that CPT violating patterns of 
neutrino spectra can be identified. It may be realized either by arrays of detectors of various 
types, or by a limited number of them with some ingenious method for analysis. Though 
highly nontrivial, this type of flavor- dependent reconstruction of supernova neutrino spectra 
is required anyway to diagnose core of the supernova in the conventional analysis assuming 
CPT invariance. The importance of the last point has been emphasized since sometime ago 
0. 

Although we formulate the problem of testing CPT with supernova neutrinos in a generic 
way, we focus on the gross (or "discrete") violation of CPT caused by different mass pat- 
terns and/or by possible large difference in mixing angle #13 in neutrinos and antineutrinos 
sectors. Then, we show that supernova neutrino can be a powerful indicator of CPT vi- 
olation. Signaling CPT violation can be done by distinguishing spectral patterns of three 
effective species of neutrinos characteristic to CPT violation, which come from unequal level 
crossing patterns at the high-density resonance of neutrinos and antineutrinos. It should 
also be noted that this type of CPT test is quite complementary to the one which measures 
small differences of neutrino and antineutrino mixing parameters. The method of looking 
for gross violation of CPT is quite insensitive to the presence of tiny difference in mixing 
parameters because the effect we are looking for is robust and only depends upon mass 
patterns. The issue of supernova neutrino as a sensitive probe for neutrino mass hierarchy 
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was first discussed in jT^. 1 



III. BASIC PROPERTIES OF SUPERNOVA NEUTRINOS AND APPROXIMA- 
TIONS INVOLVED IN THE TREATMENT 

In this section, we briefly summarize the basic properties of supernova neutrinos and 
the approximations involved in our treatment. Our description will be a very brief one 
and we refer for detailed discussions, on which our treatment and notations will be 
based. The great simplification that occurred after the work is published is that the large 
mixing angle (LMA) region of the solar Mikheyev-Smirnov-Wolfenstein (MSW) solution [2J| 
is selected out both in neutrino and antineutrino sectors by all the solar and the KamLAND 
experiments, respectively 0,0]. 

We assume three flavor mixing scheme of neutrinos with the standard form 0] of lepton 
flavor mixing matrix, the Maki-Nakagawa-Sakata (MNS) matrix 25l |. 
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where Cjj and Sjj = 1-3) imply cos and sin 6^-, respectively. The lepton mixing matrix 
U relates the flavor eigenstate to the mass eigenstate as u a = U a i^i, where a = e,fx,T and 
i = 1, 2, 3. The mass squared difference of neutrinos is defined as Amfj = mj — m| where 
rrii is the eigenvalue of the ith mass-eigenstate. To distinguish antineutrino mixing matrix 
from that of neutrinos we place a "bar" onto the corresponding mixing parameters. 

In the analysis in this paper, we restrict ourselves into the simplified ansatz for supernova 
neutrinos. That is, u^, u^, u T , and v T are treated as a single component denoted as v x . It is 
a good approximation because they interact with surrounding matter only through neutral 
current (NC) interactions, and hence they are practically physically indistinguishable with 
each other. Under the approximation, supernova neutrinos consist of the three components, 
z/ e , v e and v x . 

It is in fact very simple to compute the neutrino flux just outside supernova for a given set 
of neutrino fluxes at neutrino sphere. To do this one first draw the level crossing diagrams of 



1 The paper contains, in addition to the general statement of utility of supernova neutrinos as a tool of 
discriminating mass hierarchy, an analysis of SN1987A data which leads the authors to conclude that (in 
page 306) "if the temperature ratio r = T Wa ./Tp e is in the range 1.4-2.0 as the SN simulations indicate, 
the inverted hierarchy of neutrino masses is disfavored by the neutrino data of SN1987A unless the H 
resonance is nonadiabatic" . While it follows the spirit of the earlier analyses l^jj, EI ! our analysis using 
the three-flavor mixing framework has physics consequences quite different from the ones spelled out in 
these papers. In fact, the ansatz tested in j^] is different from the hypothesis we have tested (which 
was relatively more disfavored) due to the three-flavor treatment of the problem. We note that most of 
the criticism posed by Barger et al. [^J does not apply to our analysis because it does not rely on the 
goodness of fit in the likelihood analysis but on the credibility of the parameters obtained as a result of fit. 
The nature of this type of analysis was already made fully transparent by Jegerlehner, Neubig, and Raffclt 
|2l| in their thorough analysis done in 1996 [2^. It would be very interesting to come back to the debate 
after having supernova simulations calibrated by the high-statistics data of future galactic supernova. 



4 



neutrinos and antineutrinos, as given in Fig. ^ The characteristic feature of the supernova 
neutrino level crossing diagram, which is unique among astrophysical objects, is that there 
are two resonances, one in high and the other low density regions, corresponding respectively 
to the atmospheric and the solar Am 2 scales. (In a typical supernova progenitor, they are 
located in helium and hydrogen burning shells, respectively.) They are referred as the H and 
L level crossings in this paper. The level crossing pattern as well as their (non-) adiabaticity 
are the decisive factors of neutrino flavor conversion in supernova jifj ]. 

It should be noticed that, because we are preparing a general CPT non- invariant frame- 
work, we have to draw diagrams of v and v separately, thereby allowing the cases with 
different mass hierarchies for neutrinos and antineutrinos. Altogether there are two and 
four cases of level crossing diagram, corresponding to the normal and the inverted mass 
hierarchies in v and v sectors, respectively. Proliferation of v diagram by a factor of 2 is 
due to inability of distinguishing rh 2 > mi or rh 2 < fh\ cases. It should be noticed that one 
can adopt one of the two conventions, which are equivalent with each other: (1) rh 2 > fh\ 
and < 9 < tt/2, or (2) < 9 < tt/4 with rh 2 > fn>i oi rh 2 < %. In this paper, we take the 
latter convention. 

An enormous simplification results in the treatment of neutrino flavor transformation 
in supernova (in fact in the envelope of the progenitor star) if the two resonances, H and 
L, are approximately independent with each other. It was argued in |17] that they are, 
based on a factor of ~30 difference between Am| x and Am^ but under the assumption that 
they are identical in neutrinos and antineutrino sectors. Fortunately, thanks to the currently 
available constraints on Am 2 and Am 2 which are already rather powerful, we can argue that 
the same approximation applies even when we relax the assumption that they are identical. 
We first note that Am| x and Am 2 21 are both in the "LMA" region; They are constrained 
to be in the regions 

2 x l(r 5 eV 2 < Am 2 21 < 2 x l(r 4 eV 2 , 

l(T 5 eV 2 < |Am 2 2i| < 2 x i(T 4 eV 2 , (2) 



former by all the solar neutrino experiments pj, while the latter by the KamLAND ex- 
periment 0], both at 3 a CL. On the other hand, Am 2 ^ and Am 2 31 are constrained to 
be 

9 x 10" 4 eV 2 < | Am 2 ! | < 6 x KT 3 eV 2 , 
4.5 x 1(T 3 eV 2 < | Am 2 ! | < 2 x 10~ 2 eV 2 , (3) 



at 99% CL by the SK atmospheric neutrino data |13|. Because of the factor of about 20 
difference it can be argued quite safely that the approximation of independent H and L 
resonances applies even in our general setting which accommodates CPT violation. 

It may be appropriate to mention here that the currently available bound on possible 
CPT violation in lepton mixing angles are rather mild, as summarized in ^5[. The current 
bound on the difference between sin 2 9\ 2 for neutrino and sin 2 #12 for antineutrinos is rather 
weak [7J . Even if we assume that 9\ 2 is in the first octant, 

|sin 2 # 12 -sin 2 £ 12 | < 0.3, (4) 

at 99.73% CL. The bound obtained for 6*13 is extremely weak, | sin 2 6*13 — sin 2 6* 13 | can be 
almost unity 12l |. 
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IV. CLASSIFICATION OF PATTERNS OF SUPERNOVA NEUTRINO SPEC- 
TRA 

In this section, after recollecting compact formulas for neutrino flavor conversion in su- 
pernova (Sec. II V A|) . we present a complete classification of spectral patterns of supernova 
neutrinos (Sec. II V EJ) . The cases of reduced degeneracy by the aid of additional informations 
from accelerator and reactor experiments will be discussed in Sec. |V] The general character- 
istics of the different flux patterns will be described in Sec. EU We will discuss in Sec. IVHI 
to what extent these additional informations help to discriminate flux patterns by limiting 
the number of possibilities. 

A. Neutrino and antineutrino spectra with LMA solution 

Under the approximations spelled out in the previous section, the neutrino fluxes that 
are to reach terrestrial detectors can be written in the compact notation 0| 

F e \ ( p l-p \ 

F g = p l-p F° . (5) 

AF X j \l-p l-p 2 + p + pJ \F X / 

The above general expression holds for both the normal and the inverted mass hierarchies. 
By "normal' and "inverted" we mean Am^ > and Am^ < 0, which will be denoted 
hereafter with subscripts iV and I, respectively, Though it is known that Am^ > by 
the solar neutrino observation, there are two possible subclasses in the antineutrino sector; 
Aifijj > and Am^ < as noted in j22|. The former and the latter will be referred to 
as the reactor-normal and the reactor-inverted hierarchies, respectively. We will keep this 
distinction with use of the combined subscript as N21 and 121 (N12 and 112) for Am^ > 
(Arn^j < 0) in the case of normal and inverted hierarchies, respectively. Altogether 
there are two and four different level crossing patterns in neutrino and antineutrino sectors, 
respectively. They are depicted in Fig. ^ 

The u e and v e survival probabilities p and p's are given as 



Pn = \U el \ 2 P H P L + \U e2 \ 2 P H (l ~ Pl) + |f4 3 | 2 (l - Ph) (6) 

pN2X = \U e l\ 2 (7) 

Pni2 = \U el \ 2 P L + \U e2 \ 2 (l-P L ) (8) 

for the normal hierarchy and 

Pi = \U el \ 2 P L +\U e2 \ 2 (l-P L ) . (9) 

Pm = \U el \ 2 P H + \U e3 \ 2 (l-P H ) (10) 

Pii2 = \Usi\ 2 PhPl + \U e2 \ 2 P H (l - P L ) + l^ 3 | 2 (l - Ph) (11) 



for the inverted hierarchy. 

We note that Am^ and Am 2 2 i are both in the "LMA" region; The former is constrained 
to be in the region 2 x 10 _5 eV 2 < Am^ < 2 x 10 _4 eV 2 by all the solar neutrino experiments 
while the latter is 10 _5 eV 2 < Am 2 2 i < 2 x lCU 4 eV 2 by the KamLAND experiments, both at 
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FIG. 1: Level crossing diagrams for neutrino flavor conversion in supernova. In neutrino sector there 
are 2 diagrams which correspond to the normal (a; Am^ > 0) and the inverted (b; Anijj < 0) 
hierarchies. In antineutrino sector there are 4 diagrams which correspond to the normal (c, e; 
Am^j > 0) and the inverted (d, f; Am^ < 0) hierarchies, each doubled by two patterns of small 
mass splittings, the reactor-normal (c and d; Am^ > 0) the reactor-inverted (e and f; Am^ < 0) 
hierarchies. In CPT invariant case the diagrams (a) to (d) remain. 

3 a CL. Then, one can argue quite safely that the L level crossings are adiabatic not only in 
the neutrino but also in the antineutrino channels, Pl = Pl = 0. Then, the survival factors 
take simple forms 

Pn = \U e2 \ 2 P H + \U e3 \ 2 (l-P H ) (12) 

PN21 = \U el \ 2 (13) 
VN12 = \U e2 \ 2 (14) 
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TABLE I: The electron neutrino and antineutrino survival factors p and p are presented for adia- 
batic and nonadiabatic H level crossings. Note that the apparent duplication due to the adiabatic 
and non-adiabatic H level crossing in the columns of pi, pN2i, and pNi2 are superficial for flavor 
conversion in supernova because of no H level crossing. 

for the normal hierarchy and 

Pi = \Ue2? (15) 

PJ21 = \U el \ 2 P H + \U e3 \ 2 (l-P H ) (16) 

PI12 = \Ue2\ 2 P H +\Uez\ 2 (l-P H ) (17) 

for the inverted hierarchy. Therefore, the distinction between Am^i > and Am^ < 
cases is just interchanging \U el \ and |C/ e2 |, as expected. 

We make a short remarks on possible roles played by the earth matter effect It is well 
known that if the H resonance is adiabatic it plays no role. In the case of non-adiabatic H 
resonance, it can play a role but again it is suppressed by the factor (2Ea eai th/ Am^) sin 2 29\s, 
where a = \^2GpN e (x) is related with neutrino's index of refraction in matter with Gp and 
N e being the Fermi constant and the electron number density, respectively. An explicit 
computation reveals that the effects of the earth matter effect is small, and moreover it 
cannot lift the degeneracy between the CPT-conserving and CPT-violating cases. Therefore, 
we do not discuss it further in the present paper. 



B. General Classification of Patterns of Supernova Neutrino Spectra 

In the rest of this paper, we focus on the cases in which the H level crossing is completely 
adiabatic or non-adiabatic; We do not discuss the intermediate case in which the H level 
crossing is admixture of adiabatic and non-adiabatic transitions. By doing so we restrict 
ourselves to the two regions of #13, roughly speaking, s 2 3 > 10~ 4 (adiabatic), or s 2 3 < 10~ 6 
(non-adiabatic). We can classify the possible situation into 4 x 8 = 32 cases depending upon 

• neutrino and antineutrino mass patterns are either normal or inverted, and if Am 2 ^ > 
or Am 2 ! < in the antineutrino sector 

• the H level crossings in neutrino and antineutrino sectors are adiabatic or non- 
adiabatic, 

as given in Table HJ In each case, the neutrino fluxes F e , F g , and F x at a detector can be 
predicted for a given set of F e °, F® , and F° at neutrino sphere. 

From the viewpoint of neutrino flavor transformation, however, there are enormous de- 
generacies in the 32 cases. First of all, the duplication due to the adiabatic and non-adiabatic 
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H level crossing in the columns of pi, pN2i, and Pnu (Fig-UJa, Fig.Hk, Fig.[T^) are superficial 
because of no H level crossing. In fact, one can show from Table U that there are only 6 
different patterns of the neutrino spectra; 
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\U e3 
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\u e2 


\v = 


\Ue3 


2 


P21 


:p = 


\u e2 
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\Uel 
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P32 


:p = 


\U e3 


2 ,p= 


\Ue2 
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P33 


:p = 


\U e3 


\v = 


\Ue3 
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P22 


:p = 


\u e2 


\v = 


\Ue2 


2 



Each pattern contains several cases of v and v mass hierarchies and (non-)adiabaticity of H 
resonance. Notice, however, that all the 32 cases must be treated as different scenarios from 
particle physics point of view. Despite degeneracies in the features of flavor conversion, each 
of the degenerate scenarios sometimes has different CPT transformation properties. 

We use abbreviated notation to represent them. For example, [y : N-AD, v : I-NAD) 
implies that neutrinos have the normal mass hierarchy and the adiabatic H level crossing, 
and antineutrinos have the inverted mass hierarchy and the nonadiabatic H level crossing. 
Then, the content of each flux pattern is: 

P31 : : N-AD, v : N21-AD), (u : N-AD, v : N21-NAD), {y : N-AD, v : I21-NAD) 

P23 : {v : N-NAD, P : I21-AD), (u : N-NAD, v : 112- AD), (u : I-AD, v : I21-AD), 
{v : I-AD,P : 112- AD) , {v : I-NAD, v : I21-AD) , {u : I-NAD, v : I12-AD) 

P 21 : (u : N-NAD, v : N21-AD) , (u : N-NAD, v : N21-NAD), {v : N-NAD, v : I21-NAD) , 
(u : I-AD, v : N21-AD), (u : I-AD, v : N21-NAD), {v : I-AD, v : I21-NAD), 
(v : I-NAD, v : N21-AD) , {v : I-NAD, v : N21-NAD) , (u : I-NAD, v : I21-NAD), 

P 32 : {u : N-AD,z/ : N12-AD) , {u : N-AD, v : N12-NAD) , {v : N-AD, 9 : I12-NAD) 

P 33 : {v : N-AD, v : 121- AD), (1/ : N-AD, v : I12-AD) 

P 22 : (u : N-NAD, p : N12-AD) , (u : N-NAD, g : N12-NAD) , (y : N-NAD, v : I12-NAD) , 
(i/ : I-AD, v : N12-AD), (v : I-AD, v : N12-NAD), {v : I-AD, v : I12-NAD), 
{v : I-NAD, v : N12-AD) , {u : I-NAD, z7 : N12-NAD) , {v : I-NAD, u : I12-NAD) , (19) 

where the one with (without) underline indicates the case with (without) CPT violation. 
Several immediate comments are in order; Most notably, only the CPT violating cases are 
involved in the latter three patterns P3 2 , P33 and P 22 . Whereas the first three patterns 
-P3I) P23 and P 2 i contain CPT conserving as well as violating cases; There are only 4 CPT 
conserving cases, two in P 2 i, and one in each of P 2 3 and in P31, and the remaining 28 cases 
are CPT violating. 

Therefore, if one is able to disentangle the latter three patterns observationally, the future 
supernova neutrino detection has a potential to discover CPT violation. We will discuss this 
possibility further in the subsequent sections. In the rest of the patterns P 2i , P 23 and P31, 
with coexistence of CPT violating and CPT conserving cases, observation of supernova 
neutrinos by itself cannot signal CPT violation nor prove CPT invariance. However, there 
are possibilities that one can make stronger cases with the help of terrestrial experiments as 
we discuss in the next section. 
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V. CASE OF REDUCED DEGENERACY WITH HELP OF OTHER TYPES OF 
EXPERIMENTS 



There are several cases in which CPT violation can be signaled more easily by combining 
SN v and V observation with some other experiments. It occurs in particular in the case 
that the next generation accelerator jl^,|2^, 3(1 ^| and/or the reactor experiments are 
able to measure 6^3 , or go down to the sensitivity to establish non-adiabatic H level crossing 
[33l | . At the stage, it may be possible that some experiments can determine the neutrino 
mass hierarchy. For recent discussions, see e.g., jsij- We explore in this section what would 
be the effect of these additional inputs for uncovering CPT violation. We do not discuss 
the cases in which CPT violation is already obvious by these additional informations, for 
example, the cases such as normal neutrino and inverted antineutrino mass hierarchies, or 
the measured values of and differ with each other at a high confidence level. 



A. Detection of #13 in reactor and accelerator experiments 

If the next generation accelerator experiments and/or the reactor measurement succeed 
to detect the effect of non- vanishing 6*13 and #13, respectively, it means that Ph = Ph = 0. 
Then, the degeneracy shrinks enormously. In each pattern of masses and level crossings, the 
cases which remain are: 
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N-AD, v 


: N21-AD), 




P23 


(y 


I-AD,P: 


I21-AD), : I-AD, v 


: 112- AD), 


P21 


(y 


I-AD,P: 


N21-AD), 




P32 


[y 


N-AD, v 


: N12-AD), 




P33 


[y 


N-AD, v 


: I21-AD),(z/ : N-AD, 


v : I12-AD), 


P22 


(y 


I-AD,P: 


N12-AD). 





(20) 



Novel feature of (J2Uj) is that the patterns P 2 i and P31 now contain only the CPT violating 
and CPT conserving cases, respectively. In this case it is sufficient to exclude the patterns 
P 23 and P 3 i to establish CPT violation. 



B. No signal of #13 in future terrestrial experiments 

Suppose that, instead of positive detection which was assumed in the above, no indication 
for nonzero #13 and #13 is obtained by the next generation accelerator and reactor experi- 
ments. If it continues to be true to the extreme sensitivity reachable by neutrino factory, it 
would imply that P# = P# = 1. Then, the degeneracy again decreases enormously, but in 
a quite different way of adiabatic H resonances, 

P 21 : (u : N-NAD, v : N21-NAD), (y : N-NAD, v : I21-NAD), 
{v : I- NAD, v : N21-NAD), {v : I-NAD, v : I21-NAD), 

P 22 : (u : N-NAD, v : N12-NAD) , (u : N-NAD, v : I12-NAD) , 
{v : I-NAD, v : N12-NAD), (u : I-NAD, v : I12-NAD), 

P 23 , P31, P 32 , P 33 : no case remains (21) 
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The only two patterns, P21 and P22 are allowed. Rejection of P21 or confirmation of P22 
implies CPT violation. 

C. The normal v and v mass hierarchies 

If the neutrino and antineutrino mass hierarchies are both normal (Am^ > 0) and if the 
value of #13 is not known, only four flux patterns remain; 

P31 : (u : N-AD, v : N21-AD), {v_ : N-AD, v : N21-NAD), 

P21 : {v : N-NAD, v : N21-AD) , (v : N-NAD, v : N21-NAD), 

P 32 : {v : N-AD,z/ : N12-AD) , (u : N-AD, v : N12-NAD) , 

P 22 : (u : N-NAD, v : N12-AD), {v : N-NAD, v : N12-NAD). (22) 

Rejection of the flux patterns P31 and P21, or confirmation of the patterns P32 and P22 
establishes CPT violation. 2 

D. The inverted v and v mass hierarchies 

If the neutrino and antineutrino mass hierarchies are both inverted type (Am^ < 0) and 
if the value of #13 is not known, only three flux patterns remain. 

P 23 : {y : I-AD, v : I21-AD), [y : I-AD, v : I12-AD), 

{u : I-NAD, v : 121- AD) , {v : I-NAD, v : 112- AD) 
P21 : {v_ : I-AD, u : I21-NAD), (u : I-NAD, z7 : I21-NAD), 
P 22 : {u : I-AD, v : I12-NAD), (u : I-NAD, z7 : I12-NAD). (23) 

To single out the pattern P 2 2 appears to be the easiest way to demonstrate CPT violation. 

Under the assumptions made in Sees. IV Bl IV CI and IV Dl CPT violation, once demon- 
strated, implies the reactor- inverted antineutrino mass hierarchy, m 2 < rh\. Or, equivalently, 
#12 is in the dark side. Thus, supernova neutrinos can in principle have sensitivity to the 
#12 light-side vs. dark-side confusion. Below, we examine the cases of additional inputs 
combined. 

E. Adiabatic H resonance and the normal or the inverted mass hierarchies 

If the H resonance is adiabatic, and if the neutrino and antineutrino mass hierarchies are 
both normal only two flux patterns remain; 

P 31 : (v : N-AD, v : N21-AD), 

P 32 : {v : N-AD, v : N12-AD). (24) 



2 Notice that the logic here is even if the sign of Am^ is measured only for neutrinos, for example, we 
assume the same sign for antincutrinos and yet the analysis can signal CPT violation. 
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If the H resonance is adiabatic, and if the neutrino and antineutrino mass hierarchies are 
both inverted the allowed flux pattern is unique 

P 23 : [y : I- AD, v : 121- AD), {u : I-AD, v : 112- AD). (25) 

In this case, there is no way of telling whether CPT is violated or not in our method. 



F. Non-adiabatic H resonance and the normal or the inverted mass hierarchies 

If the H resonance is non-adiabatic, distinction between mass hierarchies, normal vs. 
inverted does not make difference in the allowed flux patterns; 

P 21 : (v : X-NAD, v : X21-NAD), 

P 22 : (u : X-NAD, v : X12-NAD) . (26) 

where X can be N (normal) or I (inverted). Note that the two flux patterns are different 
from 



VI. CHARACTERISTIC FEATURES OF PATTERN DEPENDENT NEUTRINO 
FLUXES 



We now discuss the characteristic features of neutrino spectra of e, e, and x flavors in each 
classified patterns of neutrino flavor transformation in supernova. So far we have formulated, 
in a generic way, the method for testing CPT violation with supernova neutrinos. In the 
rest of this paper, we concentrate on testing CPT violation caused by difference between 
neutrino and antineutrino mass patterns as well as their mixing angles #13 and #13. We take 
012 = 612 in the following analysis. (Note that all the mixing angles are in the first octant 
in our convention, and 823 does not come into play.) In some cases based on the Garching 
simulation we need enormous accuracies of less than a few % to distinguish between varrious 
spectral patterns. (See Sec. lVIIl ) In such cases there is no hope of establishing CPT violation 
if the effects caused by small differences in mixing angle # 12 in v and V sectors and the one 
from neutrino mass pattern coexist. 

Though the current bounds on difference between # 12 and #i 2 are rather mild ones it is 
quite possible that the room between them will be tighten up as the KamLAND experiment 
proceeds. The choice Q\2 = 612 could become mandatory if the low-energy solar neutrino 
measurement [35| and the dedicated reactor 612 experiments \[~ 



3y| are both realized. 



A. Characteristic features of spectral patterns of neutrino fluxes 

To give the readers a feeling if the neutrino spectra that arise in the six different patterns 
can be distinguished, we give an illustration using a model flux based on Livermore simula- 
tion 37]. The parameters that characterize spectral form of the flux are given in Table ITTT1 
in Sec. IVII1 where comparison between results with the other two flux models based on the 
Garching simulation is carried out. In this subsection we employ the Livermore flux because 
it is, at least, most suitable for illustrative purpose, having clear differences among spectral 
shapes of three effective neutrino species, as shown in Fig. |21 We note that the results 
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FIG. 2: The primary energy spectra of three effective neutrino species v e (red solid curve), v e 
(green dashed curve), and v x (blue dotted curve) of neutrinos just outside the neutrino sphere are 
shown. The flux model based on the Livermore simulation whose parameters are given in Table ITTTI 
in Sec. I VI II The absolute normalization is arbitrary. 

with the Livermore parameters are very similar to the ones with the pinched Fermi-Dirac 
distribution used in a vast amount of literatures, e.g., in [TtL l3sj| . 

We draw in Fig. El the spectra of the three effective neutrino species u e , i> e , and v x of six 
different patterns of P 3 i to P22 defined in (|18|). One can recognize that the six patterns of 
neutrino flavor conversion in supernova are quite different with each other partly due to the 
"optimistic" choice of the parameters. The difference in spectral patterns in three species 
of neutrinos is the key to discriminate six different scenarios of flavor transformation. For 
more quantitative understanding, we give in Table ITT1 the average values of energies (E a ) 
a = e, e and x), and the width parameter (AE a ) = a/ (E 2 ) — (E a ) 2 which is also used in 
39l ] . To make the distinction among the six patterns clearer we also give in Table |H] the 
ratios (Ep) / (E s ) and (AEp) / (AE S ) for f3 = e and x, assuming that the denominators would 
be the best determined parameters. The distinction among the different flavor conversion 
pattern is obvious. 

Notice that the upper three patterns, P31, P23, and P21 contain CPT conserving as well 
as violating cases, whereas the lower three patterns P 32 , P33, and P 2 2 consist solely of CPT 
violating ones. Therefore, the above discussion applies, upon restriction to the upper three 
cases, to the conventional CPT conserving cases as well. 

We make comments on some notable features of the results presented in Table |HI 

• CPT conserving cases 

It may be instructive to understand the feature of the CPT conserving cases contained 
in the patterns P 3 i , P23 , and P 2 i . In the pattern P 31 (the case of normal hierarchy and 
adiabatic H resonance) the H resonance is in the neutrino channel, and (E e )/{E S ) > 1 
because at the terrestrial detector is dominantly composed by u x at the neutri- 
nosphere [26(. In the pattern P23 (the case of inverted hierarchy and adiabatic H res- 
onance) the H resonance is in the antineutrino channel, and therefore (E e ) / (E s ) < 1. 
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Neutrino Energy (MeV) 



FIG. 3: The energy spectra of the three effective neutrino species v e (red solid curve), v e (green 
dashed curve), and v x (blue dotted curve) at terrestrial detectors corresponding to six different 
patterns of neutrino flavor conversion P31 to P22 defined in (|18[) are presented. The Livermore 
flux model with the same values of parameters as in Fig. |2]is used. Although an over-all absolute 
normalization is arbitrary, the relative normalization between the six patterns are meaningful. 
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P31 


24.0 


16.8 


18.6 


12.0 


8.8 


11.3 


1.43 


1.10 


1.36 


1.29 


P23 


18.7 


24.0 


18.1 


11.6 


12.0 


10.7 


0.78 


0.75 


0.96 


0.89 


P21 


18.7 


16.8 


19.6 


11.6 


8.8 


11.6 


1.11 


1.17 


1.31 


1.32 


P32 


24.0 


20.5 


17.7 


12.0 


11.1 


10.7 


1.17 


0.86 


1.07 


0.96 


P33 


24.0 


24.0 


17.1 


12.0 


12.0 


10.3 


1.0 


0.71 


1.0 


0.86 


P22 


18.7 


20.5 


17.7 


11.6 


11.2 


12.6 


0.91 


0.86 


1.04 


1.13 



TABLE II: The averaged energies and the width parameters of three species of neutrinos are 
presented by using the spectra based on the Livermore simulation with the parameters given in 
Table UTTI The upper three patterns, P31, P23, and P21 contain CPT conserving as well as violating 
cases, whereas the lower three patterns P32, P33, and P22 consist solely of CPT violating ones. 



The distinction between P31 and P23 allows to distinguish the normal and the inverted 
mass hierarchies 0|. In the pattern P 2 i (the case of non-adiabatic H resonance) v e 
(z/ e ) at the earth are superposition of v e (z/ e ) and z/j, at the neutrinosphere . These 

features are extensively discussed by many authors 0, H H H H B H E3 • 

• CPT violating cases 
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As one can recognize from Table HP there is a general tendency that (E x ) / (E s ) is small 
in CPT violating cases. Unfortunately, it does not guarantee unique identification of 
them because of the similar small ratio of P 23 . But the latter has distinctive feature 
that all the ratios (E a ) / (E s ) and (AE a ) / (AE S ) (a = e, x) are smaller than unity. The 
feature may allow unique identification of P23, and hence that of CPT violating cases 
by elimination. 

B. Approximate Analytic Expression of the Flux Composition 

To facilitate clear understanding and to complement the drawing of figures we give ap- 
proximate expressions of fluxes. They will help understanding the features of Fig. El We 
use the approximations Si 3 1 and si 3 <C 1 (assuming that the former is true) which are 
numerically valid. Notice that it does not mean that we restrict to the case of non-adiabatic 
H-level crossing. The approximation applies to the expressions of fluxes after taking Pjj = 
or Pjj = 1 etc. to merely simplify the expressions. 

• Pattern P 3i 

In this case p = |£/e 3 | 2 = sf 3 and p = \U e i\ 2 = c\ 2 c\ 3 . Then, the u e , u e , and u x spectra are 
given by: 



• Pattern P 23 

In this case p = \U e2 \ 2 = s 2 2 c i3 and p = |t/ e3 | 2 = sf 3 . Then, the u e , v e , and v x spectra are 
given by: 



4F X 



2 pO , 2 pO pO 

-2 -2 pO , (-, _ -2 -2 \ pO -2 pO , -2 pO 
c 12 c 13- r e 1 \ l L 12 c 13J r x ~ c 12 r e ^ b 12 r x> 

« + (1 - c\ 2 ^)F! + (2 + s 2 13 + 44)^° 
F e ° + s\ 2 F» + (2 + c\ 2 )Fl 



(27) 



F e 

F- e 

4F X 



S 12 C 13-^e° + (1 S 12 C 13)-^r ~ S 12-^e° + C 12-^x' 
^13-^e + C 13 F X m F x , 

(1 - s 2 12 c 2 13 )F Q e + c 2 3 F° + (2 + s\ 2 c{ 3 + s 2 13 )F° x 
c\ 2 Fl + F! + (2 + s 2 l2 )Fl 



(28) 



• Pattern P 5 
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In this case p = \U e2 \ 2 = s\ 2 c\ 3 and p = \U e i\ 2 = c\ 2 c\ 3 . Then, the u e , u e , and v x spectra 
are given by: 



F _ 2 2 p o , /-, 2 2 n f o _ 2 p o , 2 p o 

P _ ^2 -2 pO , /i _ ^2 -2 \ pO ^ ^2 pO , ^2 pO 

af x = (i - + (i - 44)^ + (2 + 4.4. + 44)^° 
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• Pattern P 32 



This case contains only CPT violating patterns. In this case p = |t/ e3 | 2 = s 2 3 and 
P = I U e2 \ 2 = sl 2 Ci 3 - Then, the u e , u e , and v x spectra are given by: 

F e = s 13 F® + c 13 F® rs F°, 

F _ _ ^2 -2 F , M _ =2 -2 ^ pO ^2 pO , ^2 pO 

4F X = + (i - s\A,)f! + (2 + 4 + 44)^° 

« + 4^° + (2 + 4)^°- (30) 



Pattern p 



This case contains only CPT violating patterns. In this case p = \U e s\ = sf 3 and 
P = \U e 3\ 2 = s 2 3 . Then, the u e , v e , and u x spectra are given by: 



Pattern P 22 



5 13-^e° + C 13-^x ~ 



F _ _ F° + r 2 F° « P° 
4P X = c 2 3 P e ° + c 2 3 F° + (2 + , 2 3 + 4)P° 

« P e ° + F° + 2F°. (31) 



This case contains only CPT violating patterns. In this case p = \U e 2\ 2 = s 2 2 c 2 3 and 
P = \U e 2\ 2 = s 2 2 c 2 3 . Then, the u e , u e , and v x spectra are given by: 

F _ „2 2 F , (-, _ 2 2 npO 2 pO , 2 pO 

F_ _ ,2 -2 pO , /i _ =2 ^ )F -2 F , -2 pO 
r e — *l2°13- r e ^ I 1 6 12 c 13j- r a; ~ ^12^6 ^ c 12- r x? 

AF X = (1 - S 2 2 c 2 3 )F e ° + (1 - 4c 2 3 )P° + (2 + S 2 2 c 2 3 + s 2 12 c 2 13 )F° x 

« c 2 2 P e ° + c 2 2 F° + (2 + s 2 2 + 4)P°. (32) 



VII. TO WHAT EXTENT CAN NEUTRINO FLUX PATTERNS BE DISCRIMI- 
NATED? 

In this section, we briefly discuss to what extent the flux patterns predicted by six cases 
from P 31 to P 22 can be discriminated observationally. Our discussion cannot be a quantitative 
one because of the lack of "standard supernova model" which has comparable accuracies 
possessed by the standard solar model. But, it may give us a feeling of how accurate should 
be the flavor-dependent reconstruction of neutrino spectra to discriminate the six different 
patterns of flavor conversion. 

We also address in Sec JVII CI the question to what extent limiting the number of possible 
flux patterns by additional inputs help identifying CPT violation. 
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A. Model dependence of the prediction to flux spectral patterns 



To reflect the best knowledges of supernova simulation currently at hand, we employ in 
this section the parametrization of fluxes which is used by the Garching group to fit their 



data HI E 



F° a (E) 



®* A 



(E a ) T(f3 a ) \{E a ) 



E 



exp -f3, 



E 



(Ea) 



(33) 



where (E a ) denotes their average energy, j3 a is a dimensionless parameter which is related to 
the width of the spectrum and typically takes on values 3.5-6. For definiteness, we assume 
(3 Ve = 3.5, (5 Vx = 4, and (3p e = 5. For the sake of comparison and to reveal dependence on 
supernova simulations we examine the three different sets of parameters, the same ones used 



m 



43[. They are the fluxes obtained from the Livermore simulation [3?]], and typical two 



model fluxes based on simulation done by the Garching group 
parameters are given in Table UTTl as L, Gl and G2. 
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The three model flux 



Model (Eo(v e )) (Eofa)) (E {u x )) 



$o(^ e ) 



L 


12 


15 


24 


2.0 


1.6 


Gl 


12 


15 


18 


0.8 


0.8 


G2 


12 


15 


15 


0.5 


0.5 



TABLE III: The parameters of the primary neutrino spectra models motivated from SN simulations 
of the Garching (Gl, G2) and the Livermore (L) groups, the same as used in [43j. We assume 
(3 Ue = 3.5, (3 Ux = 4, and /3 Pe = 5. 

In Fig. HI the spectra of the three effective neutrino species u e , v e , and v x of six different 
patterns of P 3i to P 22 defined in (jTSjl are plotted by taking model parameters Gl (upper 
figures) and G2 (lower figures) in Table 11111 As one can recognize the six flavor conversion 
patterns is much harder to distinguish than the case of Livermore flux. 

To obtain a hint on how much accuracy is needed to disentangle these six patterns we 
give in Table IIVI the averaged energies and width parameters of three species of neutrinos 
as done in Table |HJ With the Garching parameters, typically a few % accuracies for the 
determination of ratios (E a )/(E S ) and (AE a ) / (AE S ) (a = e and x) are required. Whereas 
in the case of the Livermore parameters, we may be able to do the job with the accuracies 
of -10 %. 

In addition to these quantities we have added in Table IIVI columns for the ratio of lumi- 
nosity of v e and u x to u e . (We note that even if the columns are added in Table ITT1 it is not 
informative because of the equi-partition of luminosity between three species.) It should not 
be too difficult to distinguish among the six patterns if the luminosity are measured at the 
level of a few % level (Gl) and 5-10% level (G2). 



B. How accurately can the supernova neutrino fluxes be determined? 

To give the readers some feeling we give a brief remark on how accurately determination of 
the supernova neutrino fluxes can be done in future observation of galactic supernovae. We 
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FIG. 4: The energy spectra of the three effective neutrino species v e (red solid curve), v e (green 
dashed curve), and v x (blue dotted curve) at terrestrial detectors corresponding to six different 
patterns of neutrino flavor conversion P31 to P22 defined in (|18[) are presented. The Garching flux 
model parameters Gl (upper figures) and G2 (lower figures) are used. 

must emphasize that the method for such flux reconstruction, which is of great importance 
solely from the supernova core diagnostics independent of testing CPT symmetry, is not yet 
developed to a sufficient level. It is the important problem that deserves thorough study 
to identify a minimal set of detectors which are capable of reconstructing fluxes which can 
watch supernova for a long run of at least ~50 years. 
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P23 


16.5 


18.0 


16.4 


8.8 


9.0 


8.5 


0.92 


0.91 


0.98 


0.94 


0.87 


0.83 


P21 


16.5 


16.0 


16.9 


8.8 


7.7 


8.8 


1.03 


1.06 


1.15 


1.15 


1.14 


1.17 


P32 


18.0 


17.3 


16.2 


9.0 


8.6 


8.5 


1.04 


0.94 


1.05 


0.99 


1.11 


0.91 


P33 


18.0 


18.0 


16.0 


9.0 


8.4 


8.4 


1.0 


0.89 


1.0 


0.93 


1.0 


0.80 


P22 


16.5 


17.3 


16.9 


8.8 


8.6 


8.1 


0.96 


0.98 


1.03 


0.94 


0.96 


0.95 
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P 3 1 


15.0 


15.0 


14.6 


7.5 


7.1 


7.4 


1.0 


0.97 


1.06 


1.05 


1.56 


1.27 


P23 


14.5 


15.0 


14.7 


7.4 


7.5 


7.3 


0.97 


0.98 


0.99 


0.98 


0.83 


0.77 


P21 


14.5 


15.0 


14.7 


7.4 


7.1 


7.4 


0.97 


0.98 


1.05 


1.05 


1.29 


1.33 


P32 


15.0 


15.0 


14.5 


7.5 


7.4 


7.3 


1.0 


0.97 


1.02 


0.99 


1.17 


0.89 


P33 


15.0 


15.0 


14.5 


7.5 


7.5 


7.3 


1.0 


0.97 


1.0 


0.97 


1.0 


0.73 


P22 


14.5 


15.0 


15.3 


7.4 


7.4 


5.9 


0.97 


1.02 


1.01 


0.79 


0.97 


0.94 



TABLE IV: The same as in Table ITU but with the flux models Gl (upper table) and G2 (lower 
table) based on the Garching simulation with the parameters given in Table ITTTI Two extra columns 
for luminosity ratios are added. 



Lacking such studies we restrict ourselves to accuracy expected for parameters deter- 
minable with v e observation in water Cherenkov detectors. In 3^ it was found that the 



parameters of the primary flux (E s ) and te = (E x ) / (E s ) can be determined to the accura- 
cies 1% (4%) and 1.5% (9%) at 3 a CL with Hyper-Kamiokande |!H (Super-Kamiokande), 
respectively. The accuracies found in [38] correspond to the ones of primary fluxes but we 
here assume that the similar accuracies can be expected for the terrestrial fluxes. If these ac- 
curacies can be extended to v e flux (which is, however, highly nontrivial) it may be possible 
to disentangle the flux patterns expected in the six different patterns of flavor conversion. 
We stress that though the accuracy required for flux determination for identifying CPT 
violation is quite a demanding one, they are at the same level as that required to diagnose 
interior of supernova core by means of neutrinos. 



C. How additional inputs from accelerator and reactor experiments help identi- 
fying CPT violation? 

In Sec. |V] we have discussed which flux patterns will remain if the next generation ac- 
celerator and reactor experiments measure #13 and ^13 and/or determine the v and v mass 
hierarchies. We discuss in this subsection how the reduction of numbers of possible flux 
patterns help to discriminate between CPT invariant and non-invariant cases. 

If the H level crossing is non-adiabatic, the two flux patterns P21 and P22 remain. In fact, 
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it is the general feature of the non-adiabatic H level crossing, independent of mass hierarchy, 
as we have seen in Sec. IV Bl In the Garching flux model Gl, the ratios of average neutrino 
energies and the widths of energy distributions of v e and v x to v e are higher by about 7 - 
20% in the flux pattern P 21 than in P 22 . At the same time, the luminosity ratios are higher 
in P 21 by about 20%. In the model G2, the luminosity ratios are again higher by 30-40% in 
the flux pattern P 21 than in P 22 . The average energy and width ratios do not differ much 
except for j^g-) * n which it is higher by about 30% in P 21 . 

If the H level crossing is adiabatic the results depend upon the v and v mass hierarchies. 
If they are both normal we are left with two flux patterns P31 and P3 2 (Sec. IVEjl . In the 
Garching flux model Gl, all the four ratios of average neutrino energies and their widths of 
the flux patterns P31 are higher by about 10% than those of P3 2 . The luminosity ratios also 
differ by about 20%. In the model G2, average energy and width ratios do not show any 
appreciable difference, but the luminosity ratios are higher by 30-40% in the flux pattern P31 
than in P3 2 . If the v and v mass hierarchies are both inverted there is no way of signaling 
CPT violation in our method. 

Therefore, it appears that there is a chance to uncover CPT violation within the accuracy 
which may be expected in a future detectors. Of course, we should not rely too much on 
the particular set of flux models to judge to what extent the different flux patterns are 
discriminable. But, the examination we have gone through in the above suggests that 
there are some possibilities of uncovering CPT violation by our method at least under 
circumstances helped by future terrestrial measurement. 

D. Supernova model dependence 

Here, we want to give a cautionary remark. Uncovering CPT violation along the way we 
discussed in this paper can only be claimed on the ground that the supernova simulations 
at the time of observation are reliable to certain extent. Even in the luckiest case in which 
we know the mass hierarchy and that #13 is in a region measurable by the next generation 
reactor and accelerator search, we need credibility in the flux model, roughly speaking, to 
10-20% level in the predictions to the ratios of average energy, width, and the luminosity of 
v e and u x to v e . 

Suppose that a future measurement of supernova neutrino flux strongly suggest CPT 
violation by preferring one (or more) of the flux patterns which consists only of CPT violating 
mass patterns, or by disfavoring any of the CPT invariant patterns. Then, one is tempted 
to conclude that CPT violation is signaled. The point is that the conclusion can be made 
firmer only by calibrating supernova simulation by accumulation of data gained by many 
explosions for consistency check. This point is a inherent weakness of the method of signaling 
CPT violation by supernova neutrino data. We want to emphasize, however, that it is quite 
thinkable that we will have reliable model simulation of explosion in a timely way, by which 
neutrino flux can be predicted with high accuracies. Neutrinos are the main engine of the 
explosion (they are like the pp neutrinos in the Sun) and, most probably, only the ordinary 
known physics is involved inside the core. 



20 



E. Implications to analyses with CPT invariance 



As we mentioned at the end of Sec. HI part of the formulas and the analysis in this paper 
contain some informations useful for conventional analysis with CPT, in particular in the 
context of flavor- dependent reconstruction of three fluxes. We mention only a few points 
below, leaving full exposition to possible future works. 

Look at the analytic formulas (|27jl and (J2H|) in Sec. I VI HI From these equations, we 
notice the following features: In ()27|) which applies to the case of normal mass hierarchy 
with adiabatic H resonance, the primary v e spectrum does not show up (or has a suppression 
factor of s1 3 ) in the v e and v e spectra observed at the terrestrial detectors. It means that 
to reconstruct the primary v e spectrum accurate measurement of u x spectrum together with 
those of v e and v e is mandatory. Similarly, the primary v e spectrum does not appear (or 
has a suppression factor of s^ 3 ) in the v e and v e spectra in (J25|l which applies to the case of 
inverted mass hierarchy with adiabatic H resonance. 

In these two cases it would be very difficult to reconstruct the primary u e (p e ) spectrum 
in the case of normal (inverted) mass hierarchy with adiabatic H resonance, because spectral 
measurement of v x is difficult. (See, however, fi§] for a possible way out.) We believe that it 
is one of the crucial problems in the program of flavor-dependent reconstruction of primary 
neutrino fluxes in the interior of supernova. 

The case of normal and inverted mass hierarchies with non-adiabatic H resonance, (|29|). 
may be the easiest one, relatively speaking, among the three CPT conserving patterns. It is 
because spectral measurement of u e and v e (if possible) would allow us to determine all the 
three primary fluxes if the separation of two Fermi-Dirac type distributions with different 
temperatures is possible, as illustrated in [38t | - 



VIII. CONCLUDING REMARKS 

We have discussed a method of using neutrinos from supernova to test CPT symmetry. 
Because of the possibility of having mass and mixing patterns of antineutrinos which can be 
different from neutrinos 32 different scenarios are allowed which differ in neutrino mass pat- 
terns and (non-) adiabaticity of high-density resonance. They produce six different patterns 
of supernova neutrino energy spectra at the terrestrial neutrino detectors, apart from small 
modification due to earth matter effect. Among the six patterns, three of them contain only 
the CPT violating cases. Even in the mixed cases of CPT invariance and violation addi- 
tional input on the values of #13 from reactor and accelerator experiments further enhance 
the possibility of identifying CPT violation. 

Future galactic supernovae watched by arrays of massive detectors may allow flavor- 
dependent reconstruction of three species of neutrino spectra, u e , v e and u x (a collective 
notation for 1/^, 9^, u T , and i> T ). Assuming capability of obtaining such the informations, 
which are also required to diagnose supernova core in conventional analysis with CPT, we 
have shown that one of the three CPT violating patterns may be singled out observationally. 
The help by the other measurement of lepton mixing parameters by reactor and accelerator 
experiments would help to identify CPT violating cases. Thus, we have shown that super- 
nova neutrino can be a powerful tool to detect possible gross violation of CPT symmetry 
such as different mass patterns of neutrinos and antineutrinos. We emphasize the potential 
power of the method; It may allow to disentangle different (1-3) and (1-2) mass hierarchies 
both in neutrino and antineutrino sectors in a single "bang" . 
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However, we also noticed the weakness of the method. We have just mentioned the flux 
model dependence in the previous section. Another drawback of the method concerns with 
its weakness at the precision test. In this paper, we have focused on the possibility of detect- 
ing CPT violation through identifying unequal mass patterns in neutrino and antineutrino 
sectors. If the two CPT violation in masses as well as mixing angles coexist, however, it 
would be very difficult to clearly distinguish the six different patterns of neutrino flavor 
conversion, the topics we are unable to address in this paper. Therefore, it is important 
to have stringent bound on difference between neutrino and antineutrino mixing angles in 
order for the supernova method for testing CPT to work. Similarly, the analysis would be 
very complicated if #13 is in the intermadiate region between adiabatic and non-adiabatic H 
resonance. 

Suppose that CPT violation is signaled with supernova neutrinos in the way described 
in this paper. Then, one may feel that the confirmation by using man made neutrino beam 
necessary because of the fundamental importance of CPT symmetry. Then, the question 
is; Is it possible to confirm CPT violation by the alternative methods? Fortunately, the 
answer is yes. The possibility of using two detectors at the different baseline with v beam 
only (not v) to determine the neutrino mass hierarchy is proposed some time ago [50]. 
Therefore, hypothesis of different mass hierarchies of neutrinos and antineutrinos can, in 
principle, be tested by separate measurement using the ir + and tt~ beams. Determining 
the sign of Am^ in antineutrino sector is harder to carry out, as discussed in |27j . In any 
way, distinguishing neutrino mass hierarchy is a very challenging experiment and large-scale 
apparatus is required. We emphasize, therefore, that indication of CPT violation given by 
supernova neutrinos can give a good starting point of vital search for the totally unexpected 
phenomenon of CPT violation. 
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